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ABSTRACT 


The  purpose  of  this  thesis  is  to  find  suitable  ways  to  design  feedback  compensators 
for  high  order  systems  by  using  Root-Locus  methods. 

As  a  starting  point  we  will  examine  a  motor  amplidyne  system  and  a  position  con¬ 
trol  system  that  were  previously  designed  using  Bode  methods.  Then  we  generalize  the 
method  and  extend  it  to  other  systems. 

1  he  final  subject  of  this  thesis  is  to  design  feedback  compensators  as  filters  by  using 
state  feedback  coeilicients  to  define  zeros  of  the  filter,  then  we  extend  this  idea  to  build 
cascade  filters. 
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I.  INTRODUCTION 


Modern  control  systems  may  he  compensated  by  placing  a  suitable  filter  in  either 
the  forward  (cascade)  patii  or  in  a  feedback  path. 

Mainly  feedback  systems  have  the  advantage  that  the  output  follows  the  command 
more  accurately,  so  the  steady  state  error  is  less,  the  bandwidth  is  greater,  and  the  speed 
of  response  is  faster.  Also  we  can  say  that  the  elfect  of  external  disturbances  and  internal 
changes  in  the  parameters  or  structure  of  the  system  is  significantly  reduced. 

The  conventional  and  more  common  way  to  design  a  Feedback  compensator  to 
meet  a  given  set  of  specifications  is  the  Bode  plot  method,  by  determining  a  suitable  gain 
cross-over  frequency  and  1  11  .urve. 

In  this  thesis  we  will  try  to  find  some  rules  to  do  the  same  design  by  using  Root- 
Focus  methods.  To  achieve  this,  first  we  will  examine  several  designs  then  try  to  put  the 
general  ideas  in  a  set  of  ruics. 

Simulation  studies  employed  Dynamic  Simulation  Language  (DSL)  and  Lwald 
package  in  the  IBM  3033,4381  main  frame. 
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II.  PRELIMINARY  EXAMPLES 

A.  HIGH  GAIN  MOTOR  AMPLIDYNE  SYSTEM 

To  study  the  ideas  for  feedback  compensation  first  \vc  choose  an  industrial  design 
that  r  ists  and  is  available  lor  comparison.  Figure  i  shows  the  whole  compensated 
system,  where 


Figure  I.  Original  Compensated  System 


G,(.v)  = 

n{s)  = 


A',  =  14.8 


1209.8 


s(s  4-  57.9 (>)(s  4-  0. 1 2+/F32) 

Q.  137s~Yt  +  62.5)(s  4-  71.65) 
(.s  4-  1  )(.s  4-  5.26)(s  4-  625) 


(1.) 


(2) 


(?) 


Figure  2  shows  the  open  loop  Bode  plot  for  the  uncompensated  system.  The  un¬ 
compensated  system  is  badly  unstable. 
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We  may  use  two  different  kinds  of  block  diagram  manipulation  to  analyze  this  sys¬ 
tem.  First  is  the  ordinary  manipulation  shown  at  Figure  3,  where 


Figure  3.  Ordinary  Manipulation 


K,  =  14.8 


6096553. 14(s  +  l)fs  +  5.26)(s  +  625) 

C  = - - - - - - -  (5) 

cq  s(.s  +  0.003+/).066)(5  4-  59.35)0  +  79.2)(s  +  275.45+JS45.S) 

Analysis  of  this  transfer  function  provides  the  open  loop  Bode  plot,  Root-Focus, 
magnified  Root-Locus  and  time  response  of  the  system  which  are  given  by  Figures  4  - 
7.  Since  there  is  no  pole  in  the  right  half  plane,  we  can  get  phase  margin  and  gain  margin 
from  these  plots. 
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Figure  5.  Comp  System  Root-Locus 
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REAL  AXIS 


T  he  second  type  of  manipulation  is  done  to  get  more  information  using  the  Root- 
Locus  method.  The  first  step  for  this  is  shown  at  f  igure  S.  Then  with  the  appropriate 
block  diagram  reduction  the  system  becomes  as  in  Figure  9,  where 


Figure  8.  First  Step  For  Manipulation 


Figure  9.  Final  System  For  Root-Locus  Examination 


G3(s)  = 


17905 


(s  +  579.97)(s  +  56.14)(j  -  26.9S+/45.19) 


(6) 
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0.0()92653(j  +  62.16)(s  +  71.92) 
7/1  (S+ |)(j  +  5.24j(5  +  625.33) 


(7) 

By  using  the  loop  transfer  function  Gll(s)  to  draw  the  Root-Locus  we  may  examine 
the  system  for  root  movements.  Since  the  G(s)  has  two  poles  in  the  right  half  plane,  the 
Bode  plot  cannot  be  interpreted  from  the  point  of  view  of  gain  margin  or  phase  margin, 
but  it  shows  the  stability  of  the  system. 

The  root  locus  for  the  uncompensated  system  is  shown  at  Figure  10. 

From  this  point  we  can  proceed  to  explore  the  design  procedure  of  the  designer. 
First  of  all,  since  this  design  was  done  using  Bode  design  methods,  it  will  be  helpful  to 
get  the  G1  and  I'll  magnitude  curve  which  is  shown  at  Figure  1 1  lor  the  original  system 
at  Figure  1. 

As  we  may  see  from  the  block  diagram  there  are  zeros  at  the  origin.  For  this  system 
we  need  at  least  two  zeros  to  keep  the  system  type  number  unchanged.  Also  the  spec¬ 
ifications  show  that  we  want  to  have  a  very  high  error  coefficient  for  steady  state  accu¬ 
racy  and  a  very  narrow  bandwidth.  Because  of  these  characteristics,  there  is  a  dipole  near 
the  origin.  It  is  possible  to  see  this  on  the  Bode  plot  of  Figure  1 1. 

The  designer  could  have  put  the  third  zero  very  close  to  origin  but  in  that  case 
building  the  system  might  difficult. 

Since  other  poles  are  far  away  compared  with  these  zeros  and  dipole,  these  compo¬ 
nents  will  decide  the  system  behaviour. 

Ihe  Bode  plot,  Root-Locus  and  magnified  Root-Locus  around  the  origin  for  the 
loop  transfer  function  GII(s)  are  given  by  Figures  12  -  14.  Since  we  have  roots  in  the 
right  half  plane,  the  Bode  plot  cannot  be  interpreted  for  phase  and  gain  margin  values. 

As  a  result  of  this  problem  we  may  say  that  for  high  error  coefficient  and  narrow 
bandwidth  systems  we  may  need  zeros  at  the  origin  and  a  dipole  very  close  to  the  origin. 


1! 
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r  Function  GH(s) 


B.  POSITION  CONTROL  SYSTEM 


As  a  second  example  we  chose  a  position  control  system  having  an  amplidyne  and 
a  DC  armature  controlled  motor  as  its  power  element.  Figure  15  shows  the  original 
compensated  system,  where 


Figure  15.  Original  Compensated  System 
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Figure  16  shows  the  open  loop  Bode  plot  for  the  uncompensated  system.  The  un¬ 
compensated  system  is  badly  unstable. 


Again  we  may  use  two  different  kinds  of  block  diagram  manipulation  to  analyze  this 


system. 

First  is  the  ordinary  manipulation  shown  at  Figure  17,  where 


Figure  17.  Ordinary  Manipulation 


a;  =  0.233 

12000(s2  +  6.5s  +  10) 

6  = - : - : - : - - -  (16) 
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Analysis  of  this  transfer  function  provides  the  open  loop  Bode  plot,  Root-Locus, 
magnified  Root-Locus,  and  time  response  of  the  system  which  are  given  by  Figures  18 
-  21 

Since  there  is  no  root  in  the  right  half  plane,  we  can  get  phase  and  gain  margins  from 
this  plot. 
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Figure  18.  Compensated  System  Open  Loop  Bode  Plot 
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Figure  19.  Compensated  System  Root-Locus 
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The  second  type  of  manipulation  is  done  to  get  more  information  using  the  Root- 
Locus  method.  The  first  step  for  this  is  shown  at  Figure  22. 


Figure  22.  First  Step  For  Manipulation 


Then  with  the  appropriate  block  diagram  reduction  the  system  becomes  as  in 
Figure  23,  where 


Figure  23.  Fiftal  System  For  Root-Locus  Examination 
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By  using  the  loop  transfer  function  GH(s)  to  draw  the  Root-Locus  we  may  examine 
the  system  for  root  movements.  Since  G(s)  has  two  roots  in  the  right  half  plane,  the 
Bode  plot  of  this  cascade  combination  cannot  be  interpreted  from  the  point  of  view  of 
gain  margin  or  phase  margin,  but  it  shows  the  stability  of  the  system. 

The  root  locus  for  the  uncompensated  system  is  shown  in  Figure  24.  From  this 
point  w'e  can  proceed  to  explore  the  design  procedure  of  the  designer. 

First  of  all,  since  this  design  was  done  using  Bode  design  methods,  it  will  be  helpful 
to  get  the  G  and  I  II  magnitude  curve  which  are  shown  in  Figure  25  for  the  original 
system  in  Figure  15. 

As  we  may  see  from  the  block  diagram  there  are  three  zeros  at  the  origin.  As  in  the 
first  example  we  have  to  have  at  least  one  zero  at  the  origin  to  keep  the  system  type 
number  unchanged.  Also  since  we  want  to  keep  the  error  coefficient  unchanged  we  have 
to  have  one  additional  zero  at  the  origin.  The  reason  for  the  third  zero  is  to  provide  a 
dipole  near  the  origin  to  make  the  system  have  a  ve  y  high  error  coefficient  for  steady 
state  accuracy  and  a  very  narrow  bandwidth. 

For  the  Gll(s).  the  Root-Locus  and  a  close  look  around  the  origin  to  the  Root- 
Locus  arc  given  in  Figures  26  and  27. 

Now  we  may  want  to  see  the  effect  of  removing  some  of  the  components  from  the 
compensator.  First,  if  we  remove  the  third  zero  the  Root-Locus  becomes  as  in 
Figure  2S.  The  system  is  unstable.  Now  we  may  want  to  bend  the  loci  by  moving  one 
of  the  poles  close  to  the  origin  toward  left.  If  we  start  to  move  the  pole  which  belongs 
to  the  dipole  at  s  =  -2.5;  Figure  29  shows  the  condition  when  this  pole  is  at  s  =  -50. 
The  system  is  still  unstable.  By  putting  this  pole  further  from  the  zero  at  s  =  -72.39, 
we  can  bend  the  loci  toward  the  left  half  plane.  We  may  interpret  this  as  follows :  w  hen 
we  move  the  pole  to  the  left  we  increase  the  sum  of  the  poles  and  therefore  the  sum  of 
the  roots.  This  affects  the  Root-Loci  by  moving  the  asymptote  centroid  to  left,  so  the 
loci  from  the  right  half  plane  cross  into  the  left  half  plane.  Figure  30  shows  the  Root- 
Loci  for  this  condition.  Now  we  may  want  to  examine  the  specifications  of  the  system 
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and  compare  them  with  our  specifications.  Since  the  dominant  roots  are  far  away  from 
the  origin  the  accuracy  and  damping  will  decrease,  and  the  bandwidth  of  the  system  will 
increase.  These  do  not  agree  with  our  specifications. 

After  ail  these  trials  and  sample  designs  we  may  summarize  the  results  as  a  set  of 
rules  for  design  of  feedback  compensation  : 

1.  To  keep  the  system  type  number  unchanged  put  at  least  the  same  number  of  zeros 
at  the  origin  as  the  original  system's  type  number. 

2.  To  keep  both  the  system  type  number  and  the  error  coellicient  same  unchanged, 
the  number  of  zeros  at  the  origin  should  be  N  +  1,  where  N  is  the  type  number 
of  the  uncompensated  system. 

3.  To  have  high  error  coefficient  and  narrow  bandwidth  for  steady  state  accuracy  in¬ 
clude  a  dipole  close  to  the  origin.  Put  the  zero  of  the  dipole  at  the  origin  for  ease 
in  physical  realization. 

In  the  next  chapter  we  will  apply  these  rules  to  several  different  kinds  of  problems 
and  sec  whether  they  work  or  not. 
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Figure  24.  Root-Locus  For  Uncompensated  System 
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Figure  26.  Root-Locus  For  The  Loop  Transfer  Function  GH(s) 


28 


Figure  28.  The  Same  System  With  Tuo  Zeros  at  Origin. 
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III.  DESIGNING  COMPENSATORS 

A.  NON-MINIMUM  PHASE  SYSTEMS 

Before  studying  systems  that  have  all  poles  and  zeros  in  the  left  half  plane,  we  may- 
want  to  examine  non-minimum  phase  systems.  If  we  have  a  root  or  roots  in  the  right 
half  plane  we  are  not  able  to  use  Bode  methods  for  designing  compensators,  therefore 
the  Root- Locus  method  will  be  the  only  useful  tool. 

To  illustrate  we  chose  the  following  transfer  function,  which  is  the  roll  mode  of  a 
vertical  take-oil' aircraft,  figure  31  shows  the  original  system,  where 


Figure  31.  Original  System 


G(s)  = 


77.187(.s  +  0.225  ±>0.6607) 

(5  -  0.0434)(s  +  1.4G4)(5  +  0.21  ±j0.%4A)(s  +  48.35  ±>19.34) 


(19) 


As  we  may  see  from  the  transfer  function,  there  is  a  real  pole  in  the  right  half  plane. 
Tor  this  reason  we  can  not  use  Bode  methods  to  compensate  the  system.  The  open  loop 
Bode  plot,  Root-Locus  and  magnified  Root-Locus  for  this  system  arc  given  by  Figures 
32  -  34. 

To  find  the  closed  loop  roots  we  make  the  block  diagram  manipulations.  Then  the 
closed  loop  transfer  function  becomes  : 


G,(s)  = 


77.187(s  +  0.225  ±>0.6607) 


(5  -  0.036)(j  +  0.21  ±yO.S4)(j  +  1.4741(5  +  4S.34  +/4S.72) 


(20) 
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Figure  32.  Uncompensated  System  Open  Loop  Bode 
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figure  33.  Root-Locus  For  Uncompensated  System 
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Figure  34.  Magnified  Root-Locus  For  Uncompensated  System 
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From  this  point  on  we  may  start  to  think  in  terms  of  Root-Locus  and  loop  transfer 
functions. 

Since  we  will  want  to  have  loci  from  the  pole  in  the  right  half  plane  toward  the  left, 
we  have  to  put  cither  a  zero  or  a  pole  to  the  lef  t  of  this  pole.  Putting  a  zero  at  the  origin 
will  not  help  us,  because  in  that  case  we  always  will  have  a  root  in  the  right  half  plane. 
So  we  should  put  the  compensator  poles  and  zeros  somewhere  in  the  left  half  plane. 
Choosing  the  specific  zero  and  pole  locations  is  an  engineering  decision  and  mostly  de¬ 
pends  on  the  specifications.  With  this  problem  we  will  try  to  illustrate  three  different 
combinations. 


First  put  a  zero  close  to  the  origin  and  a  pole  far  away,  both  in  the  left  half  plane. 
In  this  case  the  loop  transfer  function  becomes  as  in  Figure  35,  where 


Figure  35.  Loop  Transfer  Function  Form 


G,(s) 


77.187(s  + 0.225 +;0.6607) 

(s  -  0.036)(s  +  0.21  +  >0.S4)(s  +  1.474)(s  +  48.34  +  >48.72) 


(21) 


7/,(s) 


A,(5  +  l) 
(s  +  50) 


(22) 


The  Root-Locus  and  magnified  Root-Locus  for  the  loop  transfer  function  are  given 
in  Figures  36  and  37.  As  we  may  see  from  the  plots  there  is  a  root-locus  segment  be¬ 
tween  the  pole  in  the  right  half  plane  and  the  zero  at  -1. 

To  find  the  value  for  A',,  we  should  get  the  tabulated  points  for  the  Root-Locus  and 
see  for  which  gain  all  of  the  roots  are  in  the  left  half  plane.  For  this  particular  filter  the 
minimum  value  for  A’,  was  203.7.  We  selected  A',  =  500  arbitrarily.  For  this  value,  our 
system  becomes  as  in  Figure  38,  where 
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Figure  37.  Magnified  Root-Locus  For  Loop  Transfer  Function 
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77. 1 87(5  +  0.225  T  j0. 6607) 

°'(5)  ”  (5  -  0.0434)(5  +  1.464j(5  +  0.21  +  j0.844)(s  +  48.35  +  >19.34) 


(23) 


im 


500(5  4-  1) 
(5  +  50) 


(24) 


Figure  38.  Compensated  System 


If  we  make  the  block  diagram  manipulation,  we  get  the  system  shown  in 
Figure  39,  where 


Figure  39.  Resulting  Equivalent  G 


_ 77.187(5  +  0.225+y0.6617)(5  +  50) _ 

Jeq[S)  (5  +  0.0269)(5  +  0.235+9'0.842)(5+  1.535)(5  +  49.67)(5+48.41+./48.63)  ( 

The  Root-Locus,  magnified  Root-Locus,  open  loop  Bode  plot  and  time  response  lor 
the  compensated  system  are  given  in  Figures  40  -  43.  As  wc  may  see  from  the  Bode  plot 
the  system  is  stable  but  we  do  not  have  enough  gain. 
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Figure  42.  Bode  Plot  Without  Outside  Gain 
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Figure  43.  Time  Response  Without  Outside  Gain 
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We  may  handle  this  by  putting  another  gain  block  outside  the  minor  loop.  The  final 
system  then  becomes  as  in  figure  44. 


figure  44.  Filial  Compensated  System 


K  may  be  selected  according  to  the  specifications.  We  chose  k  =  100  lor  illus¬ 
tration. 


G(s) 


77.I87U  4-  0.225  +  j'O.6607) 

(s  -  0.0434)(.s  +  1.4(4)ls  4-  0.21  +  y'0.844H5  +  48.35  +  >49.34) 


(26) 


//,(5) 


500(5  4-  1) 
{S  +  50) 


(27) 


I  he  open  loop  Bode  plot  and  time  response  for  this  particular  values  are  given  in 
Figures  45  and  46.  We  may  change  the  loop  gain  according  to  the  specifications.  1- fleets 
of  changing  the  loop  gain  will  be  shown  for  the  third  combination. 


45 


46 


The  second  type  of  configuration  puts  two  complex  zeros  close  to  the  origin  and  two 
poles  far  away,  all  in  the  left  half  plane.  In  this  case  the  loop  transfer  function  becomes 
as  in  Figure  47,  where 


Figure  47.  Loop  Transfer  Function  Form 


G,(s)  = 


77.187(s  +  0.225  +  /).6607) 


(s  -  0.036)(s  +  0.21  +  yo.84j(s  +  1.474)(s  +  48.34  +  /4S.72) 


7/2(s)- 


A'2(s+1  +j) 

{s  +  20)(s  +  30) 


(28) 


(29) 


fhe  Root- Locus  and  magnified  Root-Locus  for  the  loop  transfer  function  are  given 
in  Figures  48  and  49.  As  we  can  see  from  the  plots  we  moved  the  loci  toward  the  left 
half  plane. 
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Figure  48.  Root-Locus  For  Loop  Transfer  Function 
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Figure  49.  Magnified  Root-Locus  For  Loop  Transfer  Function 


To  find  a  value  for  K}  ,  we  should  get  the  tabulated  points  for  the  Root-Locus  and 
see  for  which  gain  all  of  the  roots  are  in  the  left  half  plane.  For  this  particular  filter  the 
minimum  value  for  K2  was  2125.8.  We  selected  K2  =  2500  arbitrarily.  For  this  value,  our 
system  becomes  as  in  Figure  50,  where 


G{s)  = 


77.1S7(s  +  0.225  +y0.6607) 


(. s  -  0.0434)(5  +  1.464J(5  +  0.21  +;C.844)(j  +  48.35  +>19.34) 


(30) 


Figure  50.  Compensated  System 


If  we  make  the  block  diagram  manipulation,  we  get  the  system  shown  in 
Figure  51,  where 


Figure  51.  Resulting  Equivalent  G 


77.187(5  +  0.225  +/).6617)(s  +  20)(.s  +  30) 
eg  S)  ~  (5  +  0.007)(s  +  0.22+>).84)(5  +  1.427)(5  +  24.5+/'6.3)(5  +  48.87+j48.3)  (32) 
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The  Root-Locus,  magnified  Root-Locus,  open  loop  Bode  plot  and  time  response  for 
the  compensated  system  are  given  in  Figures  52  -  55.  As  we  may  see  from  the  Bode  plot 
the  system  is  stable  but  the  gain  is  too  low. 


Figure  52.  Root-Locus  For  Final  System 
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Figure  53.  Magnified  Root-Locus  For  Filial  System 
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Figure  54.  Bode  Plot  Without  Outside  Gain 
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Time  Response  Without  Outside  Gain 


We  may  handle  this  by  putting  another  gain  block  outside  the  minor  loop.  The  linal 
system  then  becomes  as  in  Figure  56. 


Figure  56.  Final  Compensated  System 


K  may  be  selected  according  to  the  specifications.  We  chose  K  =  100  lor  illus¬ 
tration. 


G(s)  = 


_ 77.IS7(.s  +  0.225  T./0.6607) _ 

(5  -  0.0434)fs  +  1.464)(s  +  0.21  +  j0.844)(s  +  48.35  +  y'49.34) 


0?) 


II2(s)  = 


2500(5  +\Tj) 

(5  +  20)(s  +  30) 


(34) 


The  open  loop  Bode  plot  and  time  response  for  these  particular  values  are  given  in 
Figures  57  and  58.  We  may  change  the  loop  gain  according  to  the  specifications.  Effects 
of  changing  the  loop  gain  will  be  shown  with  the  next  combination. 

The  third  type  of  combination  puts  two  complex  zeros  and  one  real  pole  close  to  the 
origin  and  one  pole  far  away,  all  in  the  left  half  plane.  In  this  case  the  loop  transfer 
function  becomes  as  in  Figure  59,  where 
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Figure  57.  Bode  Plot  For  Final  System  Loop  Gain  =  100 
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Figure  59.  Loop  Transfer  Function  Form 


G,(i) 


77.187(5  +  0.225  T/t.6607) 

(5  -  0.036)(s  +  U.2 1  +  y'U.S4)(5  +  1.474)(.s  +  4S.34  4^/48.72) 


(35) 


k\(s  +  2  +./0.5) 
//l<^  =  (5  +  l)(5  +  50) 


(36) 


T  he  Root-Locus  and  magnified  Root-Locus  for  the  loop  transfer  function  arc  given 
in  Figures  60  and  61.  As  we  may  see  from  the  plots  the  locus  is  moved  toward  the  left 
half  plane. 


59 


■ 

m 

■ 

u 

■ 

u 

■ 

m 

I 

■ 

* 

■ 

m 

9 

» 

m 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

-5.0  M.5  H.O  -3.5  -3.0  -2.5  -2.0  -1.5  -1.0  -0.5  0.0  0.5  I.C 

REAL  AXIS 


e  61.  Magnified  Root-Locus  For  Loop  Transfer  Function 


To  find  a  value  for  A'3,  we  should  get  the  tabulated  points  for  the  Root-Locus  and 
see  for  which  gain  all  of  the  roots  are  in  the  left  half  plane.  For  this  particular  filter  the 
minimum  value  for  K}  was  81.6.  We  selected  K 3  =  250  arbitrarily.  By  this,  our  system 
becomes  as  in  Figure  62,  where 


Figure  62,  Compensated  System 


6(s)  = 


77.187{s  +  0.225  +  />.6607) 


(5  -  0.0434)(5  4-  1.464)(s  +  0.21  T ,/0.844)(s  4-  4S.35  +  >49.34) 
250(5  4-  2  +  0.5 j) 


//3(s)  = 


(s  4-1  )(s  +  50) 


(37) 


(38) 


If  we  make  the  block  diagram  manipulation,  we  get  the  system  shown  in 
Figure  63,  where 


Gcq(s)  = 


77.187(5  +  0.225  +  /).6617)(5  4-  1)(5  4-  50) 


(5  4-  0.19)(5  4-  0.236  +  j0.87)(5  4-  1.52)(5  4-  49.84)(s  4-  48.37  +  y48.6) 


(39) 
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Figure  63.  Resulting  Equivalent  G 


The  open  loop  Bode  plot  and  time  response  for  the  compensated  system  are  given 
in  Figures  64  and  65.  As  we  can  see  From  the  Bode  plot  the  system  is  stable  but  the  gain 
is  too  low. 
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We  may  handle  this  by  putting  another  gain  block  outside  the  minor  loop.  The  final 
system  then  becomes  as  in  Figure  66. 


Figure  66.  rinal  Compensated  System 


K  may  selected  according  to  the  specifications.  We  chose  different  values  for  illus¬ 


tration. 


77.1S7(5  +  0.225  +/.I.6607) 

G('S)  =  (.s  -  0.0434)(s  +  1.464)(s  +  0.21  T >0.S44)(s  +  48.35  +  >49.34)  '  U) 


IhW  = 


25U(r  +  2  +  0.5/) 
(s  +  1  )(s  +  50) 


fhe  open  loop  Bode  plot  and  time  response  for  K=  100  are  given  in  Figures  67  and 
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Figure  67.  Bode  Plot  For  Final  System  Loop  Gain  =  100 
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We  can  change  the  loop  gain  according  to  the  specifications.  Figures  69  and  70  show 
the  time  responses  for  this  system  when  this  gain  is  500  and  1000.  As  we  may  see  from 
the  time  response  we  can  change  the  amount  of  overshoot  by  changing  this  gain. 

As  a  result  of  this  problem  we  can  summarize  the  steps  for  designing  feedback 
compensators  for  non-minimum  phase  systems  as  following  : 

1.  Find  the  closed  loop  transfer  function  and  get  the  Root-Locus  for  this. 

2.  Pick  the  filter  poles  and  zeros  according  to  the  Root- Locus  rules. 

3.  Get  the  Root-Locus  for  the  loop  transfer  function.  Tabulate  the  gain  values,  and 
select  the  gain  that  puts  all  the  roots  to  the  left  half  plane. 

4.  Rearrange  the  system.  Look  for  the  Bode  plot  and  time  response.  If  they  are  ac¬ 
ceptable,  leave  it  that  way  otherwise  select  another  gain  outside  the  minor  loop. 
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B.  FLEXIBLE  ROCKET  CONTROL  SYSTEM 


As  a  second  example  we  chose  the  control  system  of  a  flexible  rocket  which  is  shown 
in  Figure  71.  T  his  example  was  extracted  from  Ref.  1. 


Figure  71.  Original  Uncompensated  System 


(jr(.s)  is  the  transfer  function  for  the  rigid  body  and  is  defined  as  : 


Gr(s)  = 


7.21 


(s+  1.6)(s-  1.48) 


Gd(5) 


337 


(s  +  333) 


.  2750 

Gs(s)  =  t - 

s2  +  42.2s  +  2750 


(x2  +  70s  +  40O0)(s2  +  22s  +  12800) 

G^f(s)  =  r  : 

(s2  +  30 5  +  5810)(x  +  3 Ox  +  12800) 


(42) 


(43) 


(44) 


(45 


And  the  structure  is  defined  as  ; 


7 \s)  = 


0.686(s  +  53)(x  -  53)(x2  -  152.2s  4-  1450Q)(x;  +  1.53. Sx  +  14500) 

(x2  4-  x  +  605)(x2  +  45. 5x  +  2600)(x2  +  2.5  fx  +  3900)(x2  +  3.99s  +  22980) 


(46) 
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After  performing  all  calculations  the  total  system  becomes  as  in  figure  72.  Since 
the  order  of  the  numerator  and  denominator  arc  quite  large,  the  transfer  function  is  split 
into  three  parts  to  fit  on  the  page. 


Figure  72.  Combined  Uncompensated  System 

Oris)  =  00,  X  x  <7,,  where  ; 

39J34576.25(.s  4-  35  +  y'52.7)(.s  +  1 1  +yl  1 2.6) 

G'01  =  (s+  15  +  y’74.7)(.s  +  1 5  +j  1 1 2. 1  jfs  +  333)(s  +  2 1 . 1  +y4S) 

(s-  21.55  +,/21.79)(.s  +  72.76  +-;!<)Q.24)(.s-  83.1  +/>  1.82) 

G° 1  ~  is  +  22.75  +  y-i(>.29)(5  +  1.255  +  y62.44)(s  +  1.995  +yI5I.57j(.s  -  1.4SJ 

(s  +  1  1.6  +  ,24.02) 

(s  +  1.6)(s  +  0.5  +  y24.6) 

As  we  may  see  from  the  transfer  function  we  have  four  zeros  and  a  real  pole  in  the 
right  half  plane.  In  addition  to  this  most  of  the  poles  and  zeros  which  are  in  the  left  half 
plane  arc  very  close  to  the  imaginary  axis. 

i  he  Root-Locus  for  the  open  loop  uncompensated  system  is  given  by  figure  73. 

1  he  system  is  unstable.  Since  there  is  a  real  pole  in  the  right  half  plane  we  can  not  easily 
interpret  the  Bode  plot  for  this  system. 

The  next  step  is  to  find  the  loop  transfer  function  as  in  previous  examples.  By  per¬ 
forming  the  regular  block  diagram  manipulation  the  closed  loop  transfer  fur  -tion  be¬ 
comes  C/,(s)  =  (;f,(.s)  x  GJs)  x  (7,j(-vj  which  is 


(48) 

(49) 


Figure  73.  Root-Locus  For  Uncompensated  System 
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(50) 


39134576.25(5  +  35  +  >52.7)(.s  +  1 1  +,,112.6)15  +  1 1.6  +.,'24) 
Gc,W_  (5  +  333)(5  +  1.9  T;i51.5)(5+  15 +>1 12.2)(5  +  13.8 +.,74.5) 

(5  +  72.8  +./100.2)(5  -  83  T >61.S)(s  -  21.6  +  >21.8) 

~  (5  +  1.5  +  >60.4)(s  +  19.9  +  >52.6)(.s  +  1  +  >24)(.s  -  0.5  +  j4.5) 


Gf3(5)  (5 +  24.8 +>44.17) 


(52) 


From  this  point  on  we  may  start  to  think  in  terms  of  Root-Locus  design.  As  we  saw 
from  the  loop  transfer  function  we  have  a  pair  of  complex  poles  and  four  zeros  in  the 
right  half  plane.  In  addition  to  this  most  of  the  poles  in  the  left  half  plane  are  very  close 
to  imaginary  axis  and  there  are  four  zeros  in  the  right  half  plane.  So  probably  some  of 
the  roots  are  going  to  end  up  at  these  zeros.  This  will  give  us  a  gain  constraint. 

Since  we  have  poles  close  to  the  origin  and  imaginary  axis  in  the  left  half  plane, 
putting  a  zero  or  zeros  at  the  origin  would  not  help  us,  because  in  this  case  we  always 
will  have  roots  in  the  right  half  plane.  So  we  decided  to  use  a  general  tjpc  of 
compensator  which  is  two  zeros  close  to  the  origin  and  two  poles  far  away.  The  purpose 
of  the  zeros  was  to  have  loci  from  the  poles  in  the  right  half  plane  terminate  on  these 
zeros.  We  made  several  trials  with  this  type  configuration  but  none  of  them  worked. 
An  illustration  of  this  is  shown  in  Figure  74,  where  G'c(s)  is  the  closed  loop  transfer 
function  and 


Figure  74.  Illustration  of  General  Configuration 


//,( 5) 


98.6(5  +  0.46  +  >4.482) 

(5  +  1UU)(5  +  200) 


(53) 
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1MRG  AXIS 
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The  Root-Locus  for  this  case  is  given  in  Figure  75. 
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Figure  75.  Root-Locus  for  General  Type  Compensated  System 
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As  we  may  see  the  system  is  unstable,  because  roots  in  the  left  half  plane  were 
moving  faster  than  the  roots  in  the  right  half  plane. 

With  this  type  of  problem,  since  the  system  has  a  12"'  order  numerator  and  17"'  order 
denominator,  designing  the  compensator  by  just  examining  the  system  would  be  impos¬ 
sible  or  require  many  trials.  So  we  decided  to  use  the  function  minimization  subroutine 
in  the  DSL  package  in  the  mainframe.  The  parameters  used  in  the  subroutine  and  the 
function  minimization  program  are  given  in  Appendix  A  and  B. 

As  the  result  of  the  program  we  found  that  we  have  to  put  the  poles  of  the 
compensator  close  to  the  origin  and  to  the  left  of  the  zeros  which  belong  to  the 
compensator.  By  putting  these  zeros  and  poles  in  this  combination,  we  make  the  roots 
in  the  right  half  plane  move  faster.  Figure  76  shows  the  block  diagram  for  the  compen¬ 
sated  system,  where  Gc(s)  is  the  closed  loop  transfer  function  and 


//2(s)  = 


14.72l.v-b  0.875  T./2.058) 
(i  +  2+  y  S.  343) 


(54) 


The  Root-Locus  for  this  compensation  is  given  in  Figure  77.  By  getting  the  tabu¬ 
lated  root  locations,  we  see  that  for  a  Root-Locus  gain  of  0.068  we  get  all  roots  in  the 
left  half  plane.  Then  our  final  compensator  becomes  ; 


,,  ,  ,  s'  4-  1.75s  4-  5 
//2  s  =  — - 

s2  +  4s  +  73.6 


(tt) 
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The  time  response  Tor  this  case  is  given  in  Figure  78.  As  we  may  see  the  system  is 
stable  but  the  settling  time  is  a  little  long.  According  to  the  specs  we  may  want  a  faster 
svstem  than  this. 


Figure  77.  Root- Locus  for  First  FM  Result 
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Running  the  Function  Minimization  program  but  changing  £  and  wn  to  get  an  out¬ 
put  with  the  desired  settling  time,  we  got  the  coefficients  for  this  compensator  as  shown 
in  Figure  79.  Gc(s)  is  the  closed  loop  transfer  function  and 


Figure  79.  Second  compensator  Resulting  from  Function  Minimization. 


im  = 


3.11(5+  1.1531  )(.y  +  4.S78) 
(s  +  3.672  +  y'5.67) 


(56) 


The  Root-Locus  for  this  is  given  in  Figure  80.  Again  by  getting  tabulated  root  lo¬ 
cations  we  see  that  for  a  Root-Locus  gain  =  0.124  we  get  all  roots  in  the  left  half  plane. 
Then  our  final  compensator  becomes  ; 


//jW  = 


S2  +  6,031  15  +  5.625 
s2  +  7.344s  +  45.632 


(57) 


The  time  response  for  this  case  is  shown  in  Figure  81.  As  we  see  the  settling  time 
and  amount  of  overshoot  are  within  specific  limits. 
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Figure  80.  Root-Locus  for  Second  FM  Result 


Figure  81.  Time  Response  for  Second  FM  Result 
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Then  our  total  compensated  system  becomes  as  in  Figure  82,  where  G'(s)  is  the  Open 
loop  transfer  function  and 


Figure  82.  Block  Diagram  of  T  otal  Compensated  System 


//3U)  = 


(,s+  1.  i  531  )(.s-f  4,878) 
(6"  +  2.672  +  j5.67) 


(^) 


As  a  result  of  this  problem  we  may  generalize  that,  if  wc  have  components  close  to 
the  imaginary  axis  in  both  left  and  right  half  plane,  we  have  to  use  pole-zero  couples 
close  to  each  other  as  a  compensator.  In  the  next  chapter  building  filters  by  using  pole 
placement  and  state  feedback  methods  will  be  discussed. 


IV.  DESIGNING  FEEDBACK  COMPENSATORS  USING  POLE 

PLACEMENT 


To  compensate  the  systems  by  state  feedback,  we  may  use  several  computer  pro¬ 
grams  to  find  the  feedback  coefficients.  Sometimes  observing  or  feeding  back  some  of 
the  states  may  not  be  possible.  For  such  cases  we  might  have  to  build  an  estimator  to 
estimate  the  states. 

The  purpose  of  this  chapter  is  to  investigate  whether  it  is  possible  to  build  a  filter 
by  using  these  feedback  coefficients  as  zeros  and  adding  extra  poles  to  the  feedback 
path. 

To  illustrate  this  we  chose  the  following  plants. 

A.  ALL  COMPONENTS  ARE  IN  THE  LEFT  HALF  PLANE 

For  this  case  the  open  loop  transfer  function  of  the  plant  is  ; 

(J(5)  = - -  (89) 

s(s+  l)(s  +  2)(s+  10)  1  ’ 

We  want  to  design  a  feedback  compensator  (a  filter)  to  stabilize  this  system  after 
determining  the  state  feedback  coefficients.  Since  the  uncompensated  system  is  fourth 
ore.  ,  in  order  to  have  full  state  feedback  we  have  to  name  four  roots. 

In  this  thesis  we  used  the  SYS  (State  Variable  Systems)  (Ref.  2]  program. 

To  use  the  computer  program  we  have  to  determine  n-1  roots  (where  n  is  the  system 
order)  and  find  the  fourth  root  from  the  characteristic  equation  to  enter  the  program. 
Our  desired  roots  are 


•h  ,2  1  5  3  —  —4 


Then  H(s)  becomes 

//(s)  =  s3  +  6s2  +  13s  +  20 


If  we  get  the  GH(s)  to  find  the  characteristic  equation 


GII(s)  = 


1000(53  +  6s2  +  13.5  +  20) 
5J  +  1 3s1  +  32  s2  +  20s 


(60) 


(61) 


(62) 
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The  characteristic  equation  then  becomes 

s4 +  10 1 3s3  +  603 25 2  +  1 30205  +  20000  (63) 

And  the  roots  of  the  characteristic  equation  are  ; 

5,  =  — 3.9SS,  s23  =  -0.997  +  >1.996,  54  =*  -1007.02  (64) 

The  root  locations  obtained  with  this  design  are  within  0.3  %  of  the  desired  values 
and  the  4''1  root  is  at  —1007.02  . 

Now  we  may  go  to  the  SVS  program  to  find  the  state  feedback  cocllicients  which 
would  place  the  roots  exactly  at  —4,  —  1  +  >2.  —1007.02.  The  state  feedback  coefficients 
we  got  rrom  SVS  are  ; 

7/(5)  =  0.9998s3  +  5.982s2  +  12.91s  +  19.8  (65) 

These  coefficients  are  negligibly  different  from  the  ones  we  obtained  by  multiplying 
tire  desired  roots.  Tor  accuracy  we  used  the  coefficients  that  we  got  from  SVS. 

f  or  all  practical  purposes  we  may  accept  the  k}  which  is  coefficient  of  s\  as  1.  Our 
state  feedback  compensated  system  then  becomes  as  in  Tigure  83,  where  G(s)  is  the 
open  loop  transfer  function  and 


Figure  83.  State  Feedback  Compensated  System. 

//(s)  =  s3  +  5.982s2  4-  12.9 Is  +  19.8  (66) 

H(s)  in  factored  form  is 

7/(s)  =  (s  +  0.996  +  >1.993)(s  +  3.997)  (67) 

Root-Loci  for  the  loop  transfer  function  Gli(s)  are  given  in  Figure  84. 


85 


To  preserve  the  unity  feedback  we  use  block  diagram  manipulation  and  the  system 
then  becomes  as  in  Figure  85. 


Figure  85.  The  Form  of  The  System  With  Unity  Feedback. 


In  Figure  85.  G(s)  is  the  open  loop  transfer  function  and 

//,(s)  =  s(s  +  2.99+  1.993)  (68) 


To  convert  the  feedback  compensator  to  make  it  more  realizable  as  a  filter,  we  ha  .c 
to  choose  the  locations  of  the  extra  poles  quite  far  away  so  they  can  not  affect  the 
transient  behaviour  of  the  system.  If  we  put  the  extra  poles  into  the  system  form  shown 
in  Figure  S3  we  get  the  required  gain  easily.  The  new  system  then  becomes  as  in 
Figure  86,  Where  G(s)  is  the  open  loop  transfer  function  and 

K(s  +  3.99)(.s  +  0.996  T./l  .993) 

=  (s  4-  300)(s  +  400)(s  +  500) 
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Figure  86.  System  With  Extra  Poles. 


To  find  the  value  of  the  K  we  get  the  tabulated  root  locations  of  the  Root-Locus  lor 
the  loop  transfer  function  G//2(j), which  is  shown  in  Figure  87. 


O  - 


300.-950.0  -200.0  -150.0  -100.0 
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Figure  87.  Root-Locus  For  Filter  Feedback  Loop  Transfer  Function. 
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Root-locus  gain  K  to  bring  the  roots  where  we  want  them  was  A'  =  1.905459  x  10' 
.  By  using  the  same  manipulation  as  in  Figure  85  if  we  redraw  the  block  diagram  of  the 
system  to  preserve  the  unity  feedback  we  get  the  final  compensated  system  shown  in 
figure  88. 


Figure  88.  Final  Compensated  System 


19.8  x  1.905459  x  107 

(70) 

(5  +  300)($  +  400)li  +  500) 

1.905459  x  107.*(5  +  2.99  ±./1.99) 

(s  +  300)(s  +  400)(s  +  500) 

(71) 

G(s)  is  the  open  loop  transfer  function  and 

G,(.s) 

//3(s)  =  - 

I  he  time  responses  for  only  state  feedback  and  for  the  filter  compensated  case  are 
Figures  89  and  90.  As  we  may  see  from  the  transient  responses  since  all  roots  are  very 
close  to  the  desired  places  the  time  responses  are  almost  exactly  same.  The  only  disad¬ 
vantage  of  this  configuration  may  be  to  build  two  identical  sets  of  poles  to  use  in  dif¬ 
ferent  places. 

B.  NON-MINIMUM  PHASE  SYSTEMS 

As  a  second  example  we  chose  a  non-minimum  phase  system  to  check  if  the  idea 
works  for  both  cases. 

For  this  case  the  open  loop  transfer  function  of  the  plant  is 


G'(5)  - 


iooo 


j(j  —  l)(s  4-  2)(5  -)-  20) 


(72) 
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Time  Response  for  Only  State  Feedback  Compensated  System 


Figure  90.  Time  Response  for  Filter  Compensated  System 
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We  want  to  design  a  feedback  compensator  (as  a  filter)  to  stabilize  this  system  after 
determining  the  state  feedback  coefficients.  Since  the  uncompensated  system  is  fourth 
order,  in  order  to  have  full  state  feedback  we  have  to  name  four  roots. 

To  use  the  computer  program  we  have  to  determine  n-1  roots  (where  n  is  the  system 
order)  and  find  the  fourth  root  from  the  characteristic  equation  to  enter  the  program. 
Our  desired  roots  are  ; 


5|j2  —  — 1  +j  1,  S3  —  —3 


(73) 


Then  I l(s)  becomes 

Il(s)  =  s3  4-  5s2  +  8s  4-  6 

If  we  get  the  GH(s)  to  find  the  characteristic  equation 

1  QOO(s3  4-  5s2  4-  8s  +  6) 

Gil(s)  =  — - - 7 - - - 

s4  +  21s  4-  18s2  -  40s 

The  characteristic  equation  then  becomes 

s4  4-  1021s3  4-  501S.s2  4-  7960s  4-  6000 

And  the  roots  of  the  characteristic  equation  are 

S,  =  -2.967,  Sj  3  =  -0.9S8  T J1.01.  s4  =  -1016 


(74) 


(75) 


(76) 


(77) 


The  root  locations  obtained  with  this  design  are  within  1.1  %  of  the  desired  values 
and  the  4''  root  is  at  -1016. 

Now  we  may  go  to  the  SVS  program  to  find  the  state  feedback  coefficients  which 
would  place  the  roots  exactly  at  —  1  +  j  1,  —  3,  —1016. 

The  state  feedback  coefficients  we  got  from  SVS  are 

//( v)  =  0.9999s3  4-  4.995s2  +  7.952s  4-  5.956  (78) 

These  coefficients  again  are  negligibly  different  from  the  Il(s)  we  first  accepted  but 
we  used  the  ones  we  found  from  SVS. 

Tor  all  practical  purposes  we  may  accept  the  &3  which  is  coefficient  of  s5  as  1.  Our 
state  feedback  compensated  system  then  becomes  as  in  Figure  91. 
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Figure  91.  State  Feedback  to  Compensate  The  System 


G(s)  is  the  open  loop  transfer  function  and 

IJ{s)  =  s3  +  4.995s2  +  7.952s  +  5.956  (79) 

I  I(s)  in  factored  form  is 

//(s)  =  (s  +  0.995  +  j0.998)(s  +  2.998)  (80) 

The  Root-Locus  for  the  loop  transfer  function  GH(s)  is  given  in  Figure  92. 


c*t 
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To  preserve  unity  feedback  we  may  use  block  diagram  manipulation.  1  he  next  form 
of  the  system  then  becomes  as  in  Figure  93  where  G(s)  is  the  open  loop  transfer  f  unc¬ 
tion  and 


H](s)  =  s(s  +  2.494  +  1.316)  (81) 


Figure  93.  The  Form  of  The  System  With  Unity  Feedback. 


To  convert  the  feedback  compensator  to  make  it  more  realizable  as  a  filter,  we  have 
to  choose  the  locations  of  the  extra  poles  quite  far  away  so  they  can  not  affect  the 
transient  behaviour  of  the  system.  If  we  put  the  extra  poles  into  the  system  f  orm  show  n 
in  Figure  91  we  get  the  required  gain  easily.  The  new  system  then  becomes  as  in 
Figure  94. 


Figure  94.  System  With  Extra  Poles 
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G(s)  is  the  open  loop  transfer  function  and 

+  2.99S)u  +  0.995  +./Q.998) 

~  (s+  3<><>)(s  +  4U0)(.s  +  5U(.) ; 

l  o  find  the  value  of  the  K  \vc  get  the  tabulated  root  locations  of  the  Root-Locus  lor 
the  loop  transfer  function  (7//2(.s), which  is  shown  in  Figure  95. 
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Figure  95.  Root-Locus  For  Filter  Feedback  Loop  Transfer  Function. 
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The  Root-locus  gain  K  to  bring  the  roots  where  we  want  them  was 
K=  1.1257119  x  10’  .  By  using  the  same  manipulation  as  in  Figure  93  if  we  redraw  the 
block  diagram  of  the  system  to  preserve  the  unity  feedback  we  get  the  final  compensated 
system  shown  in  Figure  96. 


Figure  96.  Filial  Compensated  System 


G(s)  is  the  open  loop  transfer  function  and 

5.956  x  1.1257119  x  107 


G,(s)  = 


(s  +  300 j(s  +  400)(s  +  500) 


= 


FI 257 119  x  \()'Ms+  2.494  +./1.316) 
(s  +  300)(s  +  400)(s  +  500) 


(S3) 


(S4) 


The  time  responses  for  state  feedback  and  for  the  filter  compensated  case  are  given 
Figures  97  and  98.  As  we  may  see  from  the  transient  responses  since  all  roots  are  very 
close  to  the  desired  places  the  time  responses  are  almost  exactly  same.  The  only  disad¬ 
vantage  of  this  configuration  may  be  to  build  two  identical  sets  of  poles  to  use  in  dif¬ 
ferent  places. 

C.  EXTENDING  THE  IDEA  TO  CASCADE  COMPENSATION 

The  procedure  for  root  placement  with  the  feedback  Filter  used  only  the  loop  transfer 
function.  We  may  compensate  the  system  with  the  same  components  and  same  gain  by 
using  cascade  compensation.  If  we  investigate  the  systems,  the  characteristic  equations 
are  the  same  with  the  filter  in  either  the  feedback  path  or  in  the  forward  path.  Except 
this  time  instead  of  the  poles,  the  zeros  of  the  filter  are  the  zeros  of  the  closed  loop 
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Figure  97.  Time  Response  For  Only  State  Feedback  Compensated  System 
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transfer  function  but  the  roots  of  the  closed  loop  transfer  function  arc  in  the  same  places 
as  with  the  feedback  compensation  scheme. 

To  show  the  result  of  this  idea  we  chose  the  same  plants.  For  the  first  one  if  we  put 
the  same  filter  in  the  forward  path  the  system  becomes  as  in  Figure  99,  where 


Figure  99.  Cascade  Compensated  System 


//,(s)  = 


1000 

s(s  4-  l)(s  +  2)(s  4-  10) 

(85) 

1.905459  x  lo'(j  +  1  +  j2){s  +  4) 

(86) 

(s  +  300)(s  +  400  )(5  +  500) 

The  time  response  for  this  system  is  given  in  Figure  100.  As  we  may  see  the  time 
responses  and  oscillatory  frequencies  are  almost  identical  with  the  feedback  compen¬ 
sation  time  response  which  is  shown  in  Figure  90.  In  this  configuration  because  of  the 
derivatives  due  to  the  cascade  zeros  we  have  faster  rise  time. 

We  checked  the  roots  of  the  systems  for  both  the  feedback  and  the  cascade  config¬ 
uration.  They  were  identical  within  0.37  %  difference. 


102 


Itc 


The  second  plant  becomes  as  in  Figure  101  if  we  put  the  compensator  in  the  for¬ 
ward  path,  where 


IW)  = 


.  1000 

s(s  —  l)(s  +  2)(s  +  20) 

(87) 

1.12571 19  x  I07(s+  1  Tj)(s+  3) 

(88) 

(s+300)(s  +  400)(s  +  500) 

The  time  response  for  this  case  is  given  in  Figure  102.  Again  as  we  may  see  the  time 
response  is  essentially  the  same  as  that  of  the  feedback  configuration  which  is  show  n  in 
Figure  98.  Here  again  we  have  a  faster  rise  time  because  of  derivatives  in  the  cascade 
filter. 

As  we  may  see  from  the  Figures  99  and  100  we  put  the  filter  before  the  plant  in  both 
cases  and  we  used  the  DSL  (Dynamic  Simulation  Language)  which  behaves  the  same 
way  with  the  hardware. 

Then  we  decided  to  simulate  the  systems  with  ALCON  (Simulation  program  for 
personal  computers).  This  program  uses  arithmethic  polynomial  calculations  to  calculate 
the  step  responses  for  the  systems. 
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Figure  102.  Time  Response  For  Cascade  Compensator  (Second  example) 
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The  resulting  time  responses  from  this  program  for  both  problems  are  given  in  Fig¬ 
ures  103  and  1C4.  As  we  see  high  frequency  ripples  exist  lor  a  short  time  period  which 
we  didn't  see  in  the  DSL  simulation. 

To  investigate  the  reason  for  this  difference  we  interchanged  the  blocks  in  Figures 
99  and  100,  that  is,  we  fed  the  error  signal  into  the  G  block  and  tiie  output  from  the  G 
block  into  the  H  block.  For  this  combination  resulting  step  responses  are  given  in  Fig¬ 
ures  105  and  106  which  are  identical  with  the  time  responses  we  obtained  from  ALCON. 

For  these  DSL  simulation  cases,  in  the  first  case  the  filter  is  before  plant,  the  plant 
works  as  Low-Pass  filter  and  the  high  frequency  ripples  are  not  present,  but  in  the  sec¬ 
ond  case  the  filter  follows  plant,  and  the  ripples  appear  in  the  output. 

We  could  not  decide  which  one  of  the  simulations  is  true.  According  to  linear  the¬ 
ory,  changing  the  places  of  the  plants  should  not  affect  to  the  behaviour  of  the  system. 

To  investigate  the  reason  for  this  difference  is  left  as  a  recommendation. 

As  a  result  of  these  examples  we  may  generalize  the  design  steps  as  f  ollows  : 

1.  Select  N  -  1  desired  roots. 

2.  Multiply  them  and  find  the  H(s). 

3.  Find  GII(s)  and  get  the  characteristic  equation. 

4.  Find  the  roots  of  the  characteristic  equation.  If  they  are  acceptably  close  to  the 
desired  roots,  use  H(s)  as  a  state  feedback  directly,  otherwise  go  to  SVS  or  another 
computer  program  to  find  the  state  feedback  coefficients  to  put  the  roots  exactly 
at  the  desired  places. 

5.  Choose  the  extra  poles  as  far  away  as  possible  and  put  them  in  the  denominator 
of  Ii(s). 

6.  Draw  the  Root-Locus  for  the  loop  transfer  function  GIl(s)  (including  the  extra 
poles)  and  get  the  tabulated  data. 

7.  Select  the  appropriate  gain  to  have  the  system  roots  in  the  desired  places  (generally 
this  will  be  very  high  gain). 

8.  If  the  purpose  is  to  design  a  cascade  filter,  put  the  resulting  filter  directly  in  the 
forward  path.  Otherwise  perform  the  block  diagram  manipulation  as  shown  in  the 
examples  and  preserve  the  unity  feedback. 


r-v. 
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Figure  104.  ALCON  Simulation  of  Second  Problem. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

In  this  thesis,  the  development  of  a  procedure  to  design  the  feedback  compensators 
by  using  Root-Locus  methods  and  pole  placement  methods  has  been  presented.  Simu¬ 
lation  results  and  worked  out  examples  have  shown  that  except  for  extremely  high  order 
systems  designing  the  compensators  with  the  Root- Locus  method  can  be  used  for  any 
plant.  The  most  important  result  is  to  design  them  using  pole  placement  concepts  and 
to  be  able  to  use  them  either  in  cascade  or  in  the  forward  path. 

Simulation  results  have  also  shown  that  the  compensator  can  be  determined  by  a 
Function  Minimization  subroutine  directly  from  the  transfer  function  of  an  ideal  plant 
based  on  gi\en  specifications. 

B.  RECOMMENDATIONS 

1.  Further  research  should  investigate  methods  to  satisfy  certain  specifications  while 
designing  compensators  with  the  Root- Locus  method. 

2.  In  this  thesis,  it  has  been  shown  that  Function  Minimization  can  be  uced  to  design 
compensators,  but  the  u<er  should  write  his  own  program  according  to  the  param¬ 
eters.  A  package  program  may  be  written  for  the  main  frame  to  make  the  procedure 
interactive  and  faster. 

3.  Further  research  may  investigate  the  clFcct  of  gain  variations  on  the  root  locations 
while  using  a  pole  placement  method  so  one  can  use  directly  multiplication  of  de¬ 
sired  roots  as  H(s). 


APPENDIX  A.  CONSTRAINT  PARAMETERS  OF  FUNCTION 

MINIMIZATION 


Subroutine  HOOKE  uses  the  Hooke  and  Jeeves  "direct  search"  method  to  find  the 
local  maximum  or  minimum  of  a  multi-parameter  criterion  function,  CF.  [Ref.  3|. 

The  algorithm  evaluates  CF  at  a  base  point,  X=  (A'(l) . A'(;V)),  then  perturbs  each 

parameter  in  turn,  by  the  amount  ±  STEP(J).  and  evaluates  CF  at  each  new  point. 

If  none  of  these  points  produces  a  better  value  of  CF,  then  the  stepsizes  are  de¬ 
creased  by  a  factor  of  BETA,  and  the  process  is  repeated. 

To  use  HOOKE,  the  user  must  initialize  the  following  arguments  and.  in  the  main 
program,  CALL  HOOKE(X, STEP, N,ITMAX,CFTOL, ALPHA, BETA, CF.Q.QQ. 
W.IPRINT.MIXMAX) 

All  of  these  arguments  must  be  initialized  in  MAIN,  except  for  X.  CF,  Q,  QQ.  and 
W.  Recommended  values  for  ALPHA  and  BETA  p-e  ALPHA  =  2.,  BETA  =  0.5.  All 
of  the  arrays,  i.c..  X.  STEP,  Q.  QQ  and  \V,  must  be  declared  and  dimensioned  in  MAIN 

The  arguments  and  their  meanings  may  be  explained  as  following. 

X  =  The  array  of  N  parameter  values.  The  user  must  supply  the  initial  guesses,  ei¬ 
ther  in  the  DSL  program  or  in  MAIN. 

STEP  =  An  array  of  dimension  N  containing  the  initial  step  sizes  to  be  used  in  the 
search. 

N  =  The  number  of  parameters  (a  positive  integer,  at  most  15). 

1 1  .MAX  =  The  maximum  number  of  function  calls  to  be  performed. 

CFTOL  =  The  error  in  the  criterion  function  to  be  reached  before  the  program 
terminates  (  difference  between  the  current  value  and  the  previous  stage  value). 

ALPHA  =  The  facto:  of  (Y  -  X)  which  is  added  to  Y  to  get  XNEW;  a  number 
greater  than  or  equal  to  1. 

BETA  =  The  stepsizc  reduction  factor;  a  number  between  0  and  1. 

CF  =  The  value  of  the  criterion  function. 

Q,  QQ,  W  =  Arrays  of  dimension  N,  to  be  used  as  work  space.  They  must  be  de¬ 
clared  and  dimensioned  in  the  main  program. 

IPRINT  =  An  integer  flag.  0  for  no  intermediate  printout.  1  for  intermediate 
printout  ofX,  CF,  the  number  of  function  evaluation  and  notification  of  step-reduction. 
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M INMAX  -  An  integer  flag,  -1  searches  for 


a  minimum.  +  1  searches  for  maxi- 
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APPENDIX  B.  COMPUTER  PROGRAM  FOR  FUNCTION 

MINIMIZATION 

,vThis  program  finds  the  necessary  parameters  to  fit  the  function  to  some 
’"'special  transient  response.  By  changing  zeta  and  omgn  the  standart  time 
’"'response  may  be  generated. 

D  COMMON  /HANDJ/FLAG ,ERRFN, A1 , A2 , B1 , B2 

TITLE  Rocket  parameters 

*  FUNCTION  MINIMIZATION 

’’'Put  initial  conditions  according  to  the  previous  assuming  parameters. 
INCON  A10=2. ,A20=5. ,B10=1. ,820=100. 

’■'Define  your  arrays  to  simulate  your  system 

ARRAY  P(1),Q(2),A(1),B(2),D(5),E(5),F(1),G(3),H(9),L(11),K(2),.  .  . 

M(2),Y(3),Z(3),S(2),T(4) 

’■'Put  the  values  for  the  array  coefficients 

TABLE  A( 1 )=333.  ,B(1-2)  =  1.  ,333.  ,D(  1-5)=1.  ,92.  ,18340.  ,984000 . 

51200000.  ,E(  1-5 )  =  1.  ,60.  ,19510.  ,558300.  ,74368000.  ,. . . 

F( 1)=2750. ,G( 1-3)=1. ,42. 2,2750. ,H( 1-9)=1. ,40. 652,. . . 

2949. 22,1048754. 7,1. 86426E+08.3.  9622E+09.1.  17521E+11,.  .  . 

2.  12662E+12 , 1.  10017E+14,L(  1-11)  =  1.  ,53.  12,30506. 757,.  .  . 

1378867.  ,1. 83946E+08 ,5. 2509E+09,3. 42441E+11,2.  85168E+12,.  .  . 

1. 43994E+14 , 1. 06498E+13, -3. 42095E+14,K( 1-2)  =  13.  3,39.  9,.  .  . 

M(  1  - 2 )  =  1 .  ,30.  , P(  1  )=25.  ,Q(  1  - 2 )  =  1 .  ,100.  ,Y(1)  =  1.  ,Z(1)  =  1.  ,.  .  . 

S( 1-2)=10. ,50. ,T( 1-4)=1. ,4. ,8. ,0. 

PARAM  R=l. 

’’'Changing  zeta  and  omgn  changes  the  specifications  of  your  system 

CONST  ZETA=0. 40 ,  0MGN=6. 0 , IC=0.  0 

INITIAL 

’■'This  part  of  the  program  generates  the  standart  time  response 
TOTERR=0. 0 

C1=SQRT( 1. 0-ZETA*ZETA) 

C2=0MGN*ZETA 

C3=0MGN*C1 

PHI=ATAN(C1/ZETA) 

*  ASSIGN  INITIAL  GUESS  FOR  FUNCTION  MINIMIZATION 
’’'This  part  gives  the  starting  point  to  the  program 

IF  (FLAG.  GE. 0.  )  GO  TO  10 
A1=A1 0 

IF  (FLAG.  GE.  0.  )  GO  TO  10 
A2=A20 

IF  (FLAG.  GE.  0.  )  GO  TO  10 
B1=B10 

IF  (FLAG. GE.  0. )  GO  TO  10 
B2=B20 
10  CONTINUE 

*Put  your  variable  coefficients  here  so  program  can  change  it 
Y(2)=A1 
Y(3)=A2 
Z( 2)=B1 
Z(3)=B2 

*Write  normal  simulation  program  for  your  system 
DERIVATIVE 
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XO=TRNFR( 8,10,0.  ,H,L,X2*2. 849) 

X2=TRNFR( 0,2,0.  ,F,G,X3) 

X3=TRNFR(4,4,0.  ,D,E,X4) 

X4=TRNFR( 0,1,0.  ,A,B,X6) 

X6=15.*ERG 

ERG=ER-X9 

X9=TRNFR( 2 , 2 , 0.  ,Y,Z,X0) 

ER=R-X0 

DYNAMICS 

t  ,cStandart  time  response  is  generated  by  this  formula 

STD=1.  0-(EXP( -C2*TIME)/C1)*SIN( C3*TIME+PHI ) 

ERR=X0-STD 
ERRSQ=ERR*ERR 
TOTERR=INTGRL( 0.  ,ERRSQ) 

’•'As  a  cost  function  integral  of  error  square  is  used 
TERMINAL 

,fThis  part  finishes  the  calculation  when  min  value  is  reached 
ERRFN  =  TOTERR 
CONTRL  FINTIM=2. ,DELT=0. 01 
END 
STOP 
FORTRAN 

vFrom  now  on  just  declare  your  step  values  and  amount  of  variables 
*  MAIN  PROGRAM  FOR  FUNCTION  MINIMIZATION 

IMPLICIT  REAL*8  (A-H.O-Z) 

DIMENSION  X(4) ,STEP(4) ,Q(4) ,QQ(4) ,W(4) 

C  X(1)=K,X(2)=P  INITIAL  VALUES  IN  DSL  ROUTINE 

STEP( 1 )=0. 5D0 
STEP( 2)=5. 0D0 
STEP( 3)=0. 5D0 
STEP( 4)=10. 0D0 

«  N=4 

ITMAX  =100 
CFTOL  =  1. 0D-6 

,  ALPHA  =  2. 0D0 

BETA  =  0. 5D0 
IPRINT  =  0 
M INMAX  =  -1 

CALL  HOOKE (  X , STEP , N , ITMAX , CFTOL , ALPHA , BETA , 

*  CF,Q,QQ,W, IPRINT, MINMAX) 

STOP 

END 
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